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Through a nonperturbative analysis on a sextic triple-well potential, we reveal 
novel aspects of the dynamical property of the system in connection with Af-fold 
supersymmetry and quasi-solvability. 


1 Introduction 

Analytical understanding of nonperturbative aspects in quantum theories is a 
quite difficult problem. We will present novel features of the following quantum 
mechanical system: 

ff = y + /"(l-j¥)" + |(l-3jV). (1) 

This is a triple-well potential having three local minima at g = 0 and q ~ ±l/g 
for -C 1, see Fig. |l|. It is a simple one-dimensional quantum mechanical 
system but the dynamics and the nonperturbative analysis are highly nontriv¬ 
ial. 
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Figure 1; The form of the triple-well potential (P. 






We can find that the traditional semi-classical approach fails. In the case 
of e 7^ 0, there is a so-called bonnce solntion as the classical solntion which has 
a negative mode in the flnctnations. The negative mode contribntes non-zero 
imaginary part of the spectra in the approximation, showing instability of the 
system. Since the spectra of the model mnst be real, the instability in the 
approximation mnst be fake. In the case of e = 0, the classical solntion is 
now an instanton having no negative mode and does not canse fake instability. 
However, another difficnlty takes place when we snm np the mnlti-instanton 
contribntion with the aid of the dilnte-gas approximation. The application 
of this approximation leads to the splitting of the gronnd-state spectra dne 
to the qnantnm tnnneling, in spite of the fact that the gronnd-state of the 
potential (Q) does not degenerate in the tree-level. As we will see, both the 
difhcnlties are circnmvented with the aid of the valley method briefly explained 
in the next section. Fnrthermore, the system (|^) possesses signihcant analytic 
properties at specihc valnes of the parameter, namely, A/'-fold snpersymmetry 
and qnasi-solvability. These issnes are snmmarized in Section ^ and §. 

2 Valley method 

The main problem in qnantnm theories concerns with the evalnation of the 
Enclidean partition fnnction: Z = J j Since the evalnation cannot 

be done exactly in general, one mnst hnd ont a proper method which enables 
one to get a good estimation of the qnantity. The semi-classical approximation 
is known to be one of the most established methods. Especially, the nses 
of instantons have been sncceeded in analyzing non-pertnrbative aspects of 
varions qnantnm systems which have degenerate vacna [jl|. However, validity 
of the approximation comes into qnestion when the flnctnations aronnd the 
classical conhgnration contain a negative mode. The appearance of a negative 
mode indicates that the classical action does not give the minimnm bnt rather 
a saddle point in the fnnctional space. In this case, one may expect that 
the path-integral is dominated by the conhgnrations along the negative mode, 
which may intnitively constitnte a valley in the fnnctional space. The valley 
method is a natnral realization of this consideration. 

At hrst, we give a geometrical dehnition of the valley in the fnnctional 
space g(r) 0: 
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The above definition (j^) can be interpreted as follows; for each fixed “height” 
S[q], the valley is defined at the point where the norm of the gradient vector 
becomes extremal. Introducing an auxiliary field F{q), we can make the valley 
equation (|^) a more perspicuous form: 

jd'r' D{t, = XF{t). (3) 

where the operator D is defined by D{t,F) = S[q]/Sq{T)6q{F). It is now 
evident that any solution of the equation of motion is also a solution of the 
valley equation (||) with F{t) = 0. 

Next, we separate the integration along the valley line from the whole 
functional integration. We parametrize the valley line by a parameter R and 
denote the valley configuration by quir). Introducing Faddeev-Popov deter¬ 
minant A[(y9/{], expanding the action S[q] around the valley configuration, and 
integrating up to the second order term in the fluctuation (Pr{t) = q'('r)—gij(r), 
we finally obtain the one-loop order result: 


Z = 
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( 4 ) 


where Gr{t) is the normalized gradient vector on the valley configuration. In 
the above, det^ denotes the determinant in the functional subspace which is 
perpendicular to the gradient vector Gr{t). The valley equation ensure 
that the subspace does not contain the eigenvector of the eigenvalue A. There¬ 
fore, we can safely perform the Gaussian integrations even when we encounter 
a non-positive mode. The extension to the multi-dimensional valley, which will 
be needed when there are multiple non-positive eigenvalues, is straightforward. 


3 Valley-inst ant ons 


Let us investigate the solutions of the valley equation (^) for the system (|^). 
In the case of e = 0, the three local minima of the potential (|^ have the same 
potential value. Thus, there are (anti-)instanton solutions of the equation of 
motion which describe the quantum tunneling between the neighboring vacua: 


Qo\r-ro) = ±- 


(lo\r-To) = ±- 


( 5 ) 
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When e 7 ^ 0, the classical solutions drastically change into the so-called bounce 
solutions which cause fake instability. On the other hand, the solutions of the 
valley equation (|^) contain a continuously deformed (anti-)instanton which 
connects the two non-degenerate local minima and is called (anti-)valley- 
instanton [||. 

The solutions of the valley equation also contain a family of the con¬ 
figurations, which tends to the trivial vacuum configuration in the one limit 
and tends to well-separated valley-instanton and anti-valley-instanton config¬ 
uration in the other limit. The latter configuration is called //-valley. The 
bounce solution is also realized as an intermediate configuration of this family, 
which is consistent with the fact that the solution of the equation of motion is 
also a solution of the valley equations. For details, see the numerical result in 
Ref. 0. It turns out that the quantum fluctuation around the bounce solution 
along this family of the valley conhgurations actually corresponds the nega¬ 
tive mode. Therefore, we can separate the integration along the negative mode 
with the aid of the valley method and thus do not suffer from the problem of 
fake instability. 

There are two distinct //-(//-)valley configurations since the curvature at 
the central potential bottom (at g = 0) is different, even at the leading order 
of from the one at the side potential bottoms (at q ~ ±l/g); the //-(//- 
)valley which satisfy g(±T/2) = 0 (T ^ 1) are different from the ones which 
satisfy g(±T/2) ~ 1/g or —1/g (T S> 1). The Euclidean action of the former 
with large separation R can be calculated by the perturbative expansion in 
A ~ 0(e-2^) as follows: 

= S^^^\R) = 2S'J^^ -eR+^T - R) - + 0(6"^^), (6) 

2 g 

while the one of the latter with large separation R can be calculated in the 
same way as, 

+ ^R- e(T - R) - (7) 

2 g 

where Sq'^ denotes the Euclidean action of one (anti-)instanton Eq. (^ and 
amounts to = l/4g^. In Eqs. (0) and (0), the fourth term can be in¬ 
terpreted as the interaction term between the valley-instanton and the anti- 
valley-instanton. Therefore, the minus sign indicates that the interaction is 
attractive. 

The other type of the solutions emerges in this case, which is asymptotically 
composed of two successive valley-instantons or two successive anti-valley- 


instantons. We call them //-valley and //-valley, respectively. These conhgu- 
rations do not appear in the case of double-well potentials since they connect 
every other vacuum. The Euclidean action of them with large separation R 
can be also calculated in the same way as, 

= S^^^\R) = 2^^) + - e(T -R) + + 0{e-^^). (8) 

2 

Note that the sign of the fourth term is plus and thus the interaction between 
the (anti-)valley-instantons in this case is repulsive. It turns out that these 
interaction terms in Eqs. (|^)-(|^) play an important role in circumventing the 
breakdown of the dilute-gas approximation and in producing physically ac¬ 
ceptable results. 


4 Multi-valley-instantons calculus 


Utilizing the knowledge of the (anti-)valley-instantons and the interactions 
between them obtained previously, we can evaluate the partition function 
Z = tre“'^^ by summing over those conhgurations made of several (anti- 
) valley-instantons. The sum of the contributions from the 2n valley-instantons 
conhguration can be written as. 



where denotes the contribution of the Jacobian and the /2-independent 
part of the determinant for one valley-instanton-pair and is calculated as = 
The function Jn,nu is given by. 
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where Ri is the distance between the (22 — l)-th and 22-th (anti-)valley- 
instanton and Ri the one between the 22-th and the (22 -|- l)-th (anti-)valley- 
instanton (mod n), see Fig. 
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Figure 2: The collective coordinates Ri and Ri for a 2n valley-instantons 
configuration. 


5 Nonperturbative shifts of the spectra 

From Eqs. (|^) and (pTf), the nonperturbative contributions to the spectra are 
determined by the following equation: 



We will solve the above equation by the series expansion in a\ En^jn^ = 
^noln±‘^^i where En^ stands for the spectra corresponding to, in the limit 
g —>■ 0, the eigenfunctions of the center potential well and E^^ for the ones 
corresponding to the parity eigenstates obtained by the linear combinations of 
the eigenfunctions of the each side potential well. The results are complicated 
because the calculation needs separate treatments depending on the value of 
the parameter e. For details, see Eqs. (5.22)-(5.27) in Ref. 0. From the 
results, we hnd that the nonperturbative contribution to the spectra of the 
particular levels vanishes at the specihc values of the parameter e as follows: 

1. e={AM+ l)/3 {AT = 1, 2, 3,...): = Eg = • • • = 0 for n+ < AT. 









■ ■ ■ = 0 for ?7,_ < J\f. 


2. e= {AN- l)/3 {N = 1, 2, 3,...): E^nl = E^nl = 


6 Large-order behavior of the perturbation series 


If we evaluate the spectra E by means of the perturbation expansion in the 
coupling constant as S„o/n± = i «no/n±fi'^^ the large order 

behavior of the coefficients (r 3> 1) can be estimated via the following 
dispersion relation |Q]: 
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(r) 


1 
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From the results obtained by Eq. o. the leading contributions for r S> 1 
reads, 


®no ~ No{<^) 2'’ r + -no + ^ + ) (13) 

a<:l ~ ^„^(02'Th + 3n± + |-|€V (14) 


where A„(e)’s are some constants depending on n and e. Equations (0) and 
(IT^) show that the perturbative coefficients diverge factorially unless the pref¬ 
actor y4„(e)’s vanish. It turns out that the disappearance of the leading diver¬ 
gence takes place only when e = ±(2n -|- l)/3 (n = 1, 2, 3,...), see Eqs. (5.31) 
in Ref. H. More precisely, we obtain the following results: 


1. e= (4A/'±l)/3 (TV = 1,2,3,...): kl„^(e) = 0 for n± < N. 

2. e = -(4Ar+l)/3 (Ar=l,2,3,...): kl 2 ^o+i(e) = 0 for m^KN. 

3. e =-(47^-1)73 (AT = 1,2,3,...): kl 2 ,„o(e) =0 for mo < N. 


7 A/^-fold supersymmetry 

An A/'-fold supersymmetric quantum mechanical system of one-degree of free¬ 
dom is, roughly speaking, a system of a pair of Hamiltonians which satishes 
intertwining relations with respect to an A/'-th order linear differential operator 
P^f as follows: 


PnHj^ — — 0 , ~ ~ !*• 


(15) 




For the general treatment and discussion on A/”-fold supersymmetry, see Ref. . 

It turns out that the system (|l|) becomes A/'-fold supersymmetric when 
e = ±(4A/’± l)/3. More precisely, when e = {AM ± l)/3, the relations ( |T^ ) are 
satished for H'j^ = H and following and Pj^\ 
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(17) 


When e = —{AM ± l)/3, the relations (|T^) are satisfied for Hj^ = H and 
following and Pfj-. 
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8 Quasi-solvability 

If we dehne A/”-dimensional vector spaces by = kerP^ and = 
kerP^, we see from Eq. ([T^) that are invariant under the action of the 
A/'-fold supersymmetric Hamiltonians: C This means that both 

are quasi-solvable p. Therefore, the system (P is quasi-solvable when 
e = ±{AM ± l)/3. 

Actually, when e = {AM ± l)/3, the solvable sector of P = reads, 

= span {0+ : n = 1,..., A/”} , 0n(g) = exp • (20) 

We note that, since C L^(R), M elements of the above are the exact 
eigenfunctions of the Hamiltonian They are analytic ai Q on the 

gf^-plane and thus, (i) the perturbative expansions of them in have nonzero 
convergent radii, and (ii) there is no nonperturbative contribution on them. 









The results (i) and (ii) are also the cases for the corresponding eigenvalues, 
which explain the peculiar results for the case of e = {4J\f ± l)/3 in Sections ^ 
and namely, the disappearance of the leading divergence and the vanishment 
of the nonperturbative spectral shifts. 

Similarly, when e = —(4A/” ± l)/3, the solvable sector of if = if^ reads, 

= span {0- : n = 1,... ,7^} , (t)~{q) = exp . (21) 

In contrast to in the case of e = {4J\f ± l)/3, Vjj- (/L i^^(R-) and thus any ele¬ 
ment of the above Vjf cannot be the exact eigenfunction of the Hamiltonian 
O- However, if we expand them in power of as, 0^ (g) = , 

they are normalizable at any hnite order in and thus the elements of the 
VjO give the perturbatively well-dehned and correct eigenfunctions. Therefore, 
(iii) the perturbative expansions of them in g^ have nonzero convergent radii 
too, but (iv) there are nonperturbative contributions on them. The results 
(iii) and (iv) are also the cases for the corresponding eigenvalues, which also 
explain the peculiar results for the case of e = ± l)/3 in Section H, 

namely, the disappearance of the leading divergence and the nonvanishment 
of the nonperturbative spectral shifts. 
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